load('meshex1.mat');
x = node(:,1);
y = node(:,2);
u = sin(x).*cos(y);
showsolution(node, elem, u);